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1 The baryon budget of clusters: importance 
and scientific context 

Groups and clusters of galaxies occupy a spe- 
cial position in the hierarchy of large-scale cos- 
mic structures because they are the largest and 
the most massive (from « 10 13 M o to over 
10 15 M Q ) objects in the universe that have had 
time to undergo gravitational collapse. The 
large masses of clusters imply that their con- 
tents have been accreted from regions of « 8 - 
40 comoving Mpc in size and should thus be 
representative of the mean matter content of the 
universe [1J. Thus, in contrast to galaxies, clus- 
ters are expected to retain the cosmic fraction of 
baryons within radii accessible to current and 
next generation observations in sub-mm, opti- 
cal/NIR and X-ray wavelengths, even if their 
diffuse gas was significantly heated by super- 
novae (SNe) and active galactic nuclei (AGN). 

The gaseous atmospheres of clusters (the in- 
tracluster medium or ICM) therefore contain a 
wealth of information about galaxy formation 
processes, including the efficiency with which 
intergalactic gas is converted into stars and the 
chemical and thermodynamic al effects on the 
evolution of galaxies and the surrounding inter- 
galactic medium (IGM) driven by ensuing feed- 
back mechanisms. Suppose, for example, that 
AGN feedback significantly heats the ICM in 
the central regions and the gas expands. The 
density of the gas in the center would be lower, 
but we should still observe the expanded gas in 
the outer regions of the cluster and should be 
able to observe the signature of heating in its 
entropy El OH. Likewise, if a certain fraction 
of ICM gas has cooled in the past, we should 
be able to observe these baryons either in the 
form of cold gas or in the form of stars 0. 
Some fraction of gas processed by galaxies is 
returned to the intracluster medium via winds 
and ram pressure stripping, as is evidenced by 
the high abundance of heavy elements in the in- 
tracluster plasma. The distribution of metals as 
a function of radius provides a record of the 



baryons recycled through galaxies and the ex- 
change of material between galaxies and sur- 
rounding gas. Multi- wavelength, holistic stud- 
ies of baryon contents of groups and clusters as 
a function of radius thus shed light on many key 
processes accompanying galaxy formation, one 
of the central unsolved problems of modern as- 
trophysics. Most of these processes are difficult 
to observe or constrain by other means. 

Knowledge of the baryon content of clus- 
ters is also a key ingredient in the use of clus- 
ters as cosmological probes. First, it is the ob- 
served stellar and hot gas components of clus- 
ters that are most often used as proxies for their 
total mass, which can be connected to theo- 
retical predictions and used to constrain cos- 
mology via evolution of the cluster abundance 
and spatial distribution 0. Second, measuring 
the fraction of the total mass in baryons within 
clusters is one of the most powerful methods 
of measuring the mean matter density of the 
universe [Q3 13 and is important in the interpre- 
tation of the cluster signal contribution to the 
small-scale CMB anisotropics in terms of the 
power spectrum normalization (H |9l OH flTTl . 
Third, since the mass fraction of hot gas mea- 
sured from observations depends on the lu- 
minosity distance at the cluster redshift, trac- 
ing evolution of this fraction as a function of 
redshift is a powerful independent geometrical 
constraint on the dark energy content of the uni- 
verse O ECU M EES M El EQ. ah of the 

above measurements require good understand- 
ing of the evolution of the stellar and hot gas 
fractions in clusters. 

Study of the baryon budget of groups and 
clusters across cosmic time is thus a truly 
important scientific theme with a wealth of 
connections to other areas of astronomical 
research, such as cosmology, galaxy formation 
and evolution, AGN physics, etc., and is 
therefore of tremendous scientific interest. 
During the past two decades, the advent of 
new observational facilities and theoretical 
models have allowed us to uncover some very 



interesting puzzles in this area (see § [2]) and 
to formulate specific questions, which we are 
poised to address in the next decade: 

> Does the baryon fraction of groups and 
clusters reflect that of the universe as a 
whole at radii and epochs probed by ob- 
servations ? 

> What are the radial distributions of differ- 
ent mass components of clusters ( stars and 
cold gas, diffuse hot gas, and dark mat- 
ter)? 

> How does the observed baryon budget in 
clusters and groups constrain models of 
galaxy formation and feedback? 

The radial distributions of all main mass 
components (stars, cold, warm, and hot gas and 
total mass) at different redshifts out to the virial 
radius in the same clusters are the main re- 
quired observational measurements. During the 
next decade sensitive multi- wavelength obser- 
vations, capable of accurately probing all the 
major matter components of clusters, should al- 
low us to tackle or make substantial progress on 
all of the outlined questions. Increasingly so- 
phisticated simulations of clusters will explore 
a range of plausible scenarios and assump- 
tions about the physics of galaxy formation and 
ICM and their effect on resulting distribution 
of baryons and total mass. Multi- wavelength, 
comparative studies of real and simulated clus- 
ter samples would allow us to use clusters as 
veritable astrophysical laboratories for study- 
ing galaxy formation and testing our theoretical 
models of structure formation and underlying 
assumptions about fundamental physics gov- 
erning the universe^ At the same time, reliable 

1 The Bullet cluster, the system in which hot gas in one of 
the merging clusters is clearly offset from center of total 
mass and galaxies, is a fine example of a system where 
multi-wavelength mapping of different components of 
cluster region provides powerful information about un- 
derlying physics of the system and properties of dark 
matter [19 1 . 



and detailed knowledge of the baryon content 
of clusters would provide the necessary basis 
for their use as precision cosmological probes. 

2 Current status 

2.1 Observations. Over the past decade, dra- 
matic progress has been made in observations 
of all the main components of clusters. Deep 
high-resolution wide-field imaging has dramat- 
ically improved the accuracy of strong and 
weak lensing measurements of cluster mass 
profiles and the number of systems with lensing 
data has been increasing dramatically [|20ll2T| . 
Such measurements provide a powerful inde- 
pendent method of measuring the distribution 
of total mass in clusters. Statistical measure- 
ments of the shear profiles around large num- 
bers of clusters in the SDSS survey have deter- 
mined the average mass profiles around clus- 
ters of different richness with exquisite preci- 
sion out to ~ 10 Mpc 11221 . 

Deep imaging in multiple optical bands, in 
particular in the near infrared (NIR), has per- 
mitted a systematic census of the stellar compo- 
nent in clusters [|23ll24ll . including the low sur- 
face brightness intracluster light (ICL). At the 
same time, measurements of stellar kinematics 
with the panoramic integral field spectrograph 
SAURON |25l have provided the first direct 
constraints on the mass-to-light ratios of stellar 
populations of some nearby cluster galaxies. 

New, high-resolution, sensitive X-ray obser- 
vations by Chandra and XMM-Newton have 
revolutionized our knowledge about properties 
of the hot ICM gas by reliably mapping gas 
density, temperature, and metallicity out to a 
significant fraction of virial radius ll26l 1271 . 
These observations show that the ICM is not 
isothermal and that the widely used /?-model 
does not describe the radial distribution of the 
hot gas density in the outer regions of clusters. 
Qualitatively better data about the thermal state 
and cooling of the gas, especially in cluster 
cores, has revealed both new insights and puz- 
zles, such as lack of very cold gas in the centers 



Il28l and high gas entropy in the outskirts 11291 . 
Deep, targeted X-ray observations of carefully 
selected relaxed clusters have determined the 
total mass profiles of groups and clusters under 
the assumption of hydrostatic equilibrium and 
put constraints on the halo concentrations and 
cuspy inner density profiles predicted by Cold 
Dark Matter (CDM) structure formation simu- 
lations 

Measurements of the Sunyaev-Zeldovich ef- 
fect (SZE) with interferometric techniques have 
produced independent constraints on the ther- 
mal state of the gas in dozens of clusters lfT6l 
[321 . At the same time, technological break- 
throughs in building large bolometric detector 
arrays have improved the mapping speed of 
the SZE receivers by an order of magnitude 
1133113411 . Such instruments now provide unique 
constraints on the ICM gas distribution in the 
outermost regions of clusters, to the virial ra- 
dius and beyond. 

Deep 21 cm imaging of the Virgo cluster 
has constrained the amount of cold HI gas 
ll35ll . while optical and molecular observations 
have revealed the presence of some cold gas in 
cluster centers (351 EU [381 ED- These mea- 
surements have started to elucidate the way in 
which the interstellar medium of galaxies gets 
stripped by the hot ICM as they enter clus- 
ters and the balance of heating and cooling 
in cluster cores. UV spectroscopy of back- 
ground quasars has given first constraints on 
the amount and properties of warm gas (T ~ 
few x 10 5 K). 

Figure [TJ shows three representative exam- 
ples out of many exciting developments out- 
lined above: measurement of the contribution 
of the central galaxy and ICL stars to the to- 
tal stellar budget of clusters out to the virial ra- 
dius, recent high signal-to-noise measurements 
of the SZE signal out to the virial radius, and 
measurements of ICM gas density profiles and 
total mass profiles from deep targeted Chandra 
X-ray observations of relaxed clusters. The fig- 
ure demonstrates that current observations are 
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Fig. 1 — Distribution of mass fraction of stars in the 
brightest cluster galaxies and ICL in clusters of different 
mass (top, 1 24 1), accurate SZE profile measurement for 
cluster AS 1063 out to the virial radius, m 2rsoo (middle, 
SPT team, in prep.), and hot gas and total mass profiles 
derived from deep Chandra data for a sample of relaxed 
clusters (bottom, ||26l ). 



becoming capable of mapping profiles of differ- 
ent cluster components to a significant fraction 
of the virial radius. 

2.2 Theory. During the last two decades self- 
consistent cosmological simulations of cluster 
formation following the dynamical evolution of 
collisionless dark matter and diffuse baryons 
in the non-radiative regime and including var- 
ious processes accompanying galaxy formation 
have improved in mass resolution and dynamic 
range by about 4 and 2 orders of magnitude, 
respectively. Despite the fact that many as- 
pects of galaxy formation remain poorly un- 
derstood, notable successes of the simulations 
include correct predictions of the gas density 
and temperature profiles outside of cluster cores 
001 SB 03, and the radial distribution of clus- 
ter galaxies [|43ll . Prediction of these quantities 
in an ab initio model of ACDM structure for- 
mation is convincing evidence that modeling of 
cluster formation is on a firm theoretical foot- 
ing. Cosmological simulations of clusters have 
demonstrated that the baryon fraction outside 
cluster cores is expected to be within « 10% 
of the universal baryon fraction Q, b . [|4H l5ll4*4l. 
However, stellar and hot gas fraction profiles do 
depend sensitively on the details of the physical 
processes included in simulations. ll2ll5ll4"5ll3Tl. 

Figure [2] illustrates the current observational 
measurements of stellar and hot gas mass frac- 
tions, and their sum - the total baryon frac- 
tion, within r5oo (the radius enclosing an over- 
density of 500p c , p c = 3H(z) 2 /SnG, r 500 ~ 
0.5r v ; r ), and compares these measurements to 
high-resolution simulations of clusters with dis- 
sipative physics of galaxy formation [|42l . 

Stellar mass fractions in observed groups and 
clusters (bottom panel) steadily decrease with 
increasing cluster mass. Clusters from cosmo- 
logical simulations shown in the figure exhibit 
at best a very weak trend. The actual values 
of /* for groups (M < 10 14 M o ) are in rea- 
sonably good agreement. For massive clusters, 
however, the /* in simulated clusters is a fac- 
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Fig. 2 — The stellar (bottom), gas (middle) and to- 
tal baryon (top) mass fractions in units of the mean 
baryon fraction (Qb/O m ) in cosmological cluster sim- 
ulations with (blue) and without (red) galaxy forma- 
tion (points with no error bars) and in observed clusters 
E6ll24ll46ll47l (points with error bars), within the radius 
r5oo enclosing overdensity of 500 times the critical den- 
sity. The figure shows fairly good agreement of observed 
and predicted gas fractions, but large discrepancies in the 
stellar fractions. The total baryon fractions of most of the 
observed clusters are significantly below the theoretical 
expectation. 

tor of ~ 2-5 larger than observed. This is 
an illustration of the well-known "overcooling 
problem .'[*] 

Surprisingly, despite the glaring discrepancy 
in /*, the comparison of hot gas fractions in the 
middle panel shows rather good agreement be- 
tween simulations and data, both in the average 
trend and in the actual values. The top panel 
shows comparisons of the total baryon fractions 



This problem should be distinguished from the "cool- 
ing flow" problem E51 . The latter refers to the lack of 
cold gas in cluster centers, while the overcooling prob- 
lem refers to a too high fraction of cooled, condensed gas 
produced in simulations over the entire formation history 
of clusters, including all of its progenitors. 



in simulations performed both in non-radiative 
regime (red triangles) and with galaxy forma- 
tion^] (blue circles) to observations. The total 
baryon fractions of both observed and simu- 
lated clusters do not show a detectable trend 
with mass. For observed systems this, by it- 
self, is a striking fact as gas in the shallower 
potential wells of groups is expected to be more 
susceptible to AGN feedback. Another note- 
worthy fact is that observed baryon fractions 
are lower than the universal value and the val- 
ues predicted by simulations for this radius 
HUBS EH. Although there is a scatter in val- 
ues of /bar from cluster to cluster, a large frac- 
tion of them is significantly below the universal 
value measured by the WMAP satellite. 

There are several different potential explana- 
tions for the low observed baryon fractions, al- 
though none of them have been convincingly 
proven yet. One key theoretical issue is system- 
atic errors in stellar masses derived from stellar 
population synthesis models I15T1I521 . but other 
unknown systematics in estimates of stellar and 
gas masses are possible. The X-ray estimates 
of the latter, for example, depend on the Hubble 
constant as / gas oc h~ L5 so a value of h = 0.65 
(2cr lower than h = 0.72 - the best fit value to 
the WMAP 5 yr data) would explain a signifi- 
cant fraction of the difference. 

The discussed comparisons of expectations 
and observations show that our current "vision" 
and understanding of the baryon budget within 
virial radius are not yet as clear as we would 
like them to be. The possibility remains that 
we may be missing key details of the picture. 
However, the situation can be dramatically im- 
proved in the next decade with the advent of a 
new generation of surveys and instruments at 
different wavelengths. 



The /b^ values in these two types of simulations are 
close to the universal but are somewhat different. Phys- 
ical processes behind these differences are understood, 
but their explanation is beyond the scope of this paper. 



3 Studying baryon budget in the next decade 

3.1 Observations. The next generation of ob- 
servations to constrain baryon budget of groups 
and clusters should cover a representative range 
of redshifts (0 < z < 2). This would allow to 
study formation of these systems as it occurs 
(and where sensitivity to dark energy is great- 
est in cosmological tests). Additionally, obser- 
vations should probe a wide range of masses, 
as smaller systems merge and are incorporated 
into larger ones, so their properties are key 
to understanding of the baryon content of the 
larger systems. Lastly, it is very important to 
trace mass components as far in radius as pos- 
sible, as cluster outskirts may contain key in- 
formation both about accretion of matter onto 
clusters and about the prior thermal history of 
baryons. The observational program can ad- 
dress the scientific questions formulated in § [T] 
by a combination of statistical studies of clus- 
ters in large wide-area surveys to study mean 
trends and deep targeted observations of indi- 
vidual systems to study variations as a function 
of specific physical properties. 

Stellar content. To obtain robust stellar mass 
measurements for cluster galaxies and ICL, 
deep, multi-band photometry is needed, with 
filters spanning the peak in the stellar compo- 
nent of the SED. This implies a combination 
of optical and infrared imaging for the redshifts 
< z < 2. Such imaging can be carried 
out efficiently using a combination of a large 
wide area optical ground-based or space imag- 
ing survey, and NIR imaging with JWST from 
space. Large area multi-band optical surveys 
could constrain ICL in low-z clusters for which 
X-ray and SZ data are available by stacking 
clusters within a given range of total mass. In- 
dependent constraints on the mass-to-light ra- 
tios of a large number of cluster galaxies with 
IFU spectroscopy on large telescopes would 
usefully extend upon the work of SAURON 
project in the recent years. Development of 
better SPS models and stellar mass estimators 



and good understanding of the associated sys- 
tematic error is a must (52]. Spectroscopic 
measurements of velocities in cluster fields in 
conjunction with theoretically-motivated mock 
catalogs could constrain the 3d distribution of 
cluster galaxies. Efficient execution of such 
a spectroscopic measurement program requires 
a highly multiplexed spectrograph capable of 
dense spatial sampling (e.g. an instrument like 
GISMO+IMACS at Magellan) or a wide-area 
IFU, on an 8 - 10 m or larger telescope. 

Hot ICM gas. Targeted deep observations 
of individual groups and clusters by Chandra, 
XMM-Newton, and Suzaku will continue to 
provide measurements of gas density and tem- 
perature profiles out to r 500 - r vir . The num- 
ber of possible targets will increase dramati- 
cally after completion of the currently ongoing 
X-ray and SZ surveys and those planned for the 
near future. It would be extremely helpful if 
by the end of the next decade a more sensi- 
tive, low background, high-resolution (< 10 ") 
X-ray telescope capable of mapping the gas dis- 
tribution to larger radii (to r v ; r and perhaps be- 
yond) at < z < 2, such as 1X0 (see white 
paper by Vikhlinin et al.), would become avail- 
able. Sensitive X-ray observations would help 
to trace evolution of hot gas fraction very accu- 
rately over the main epochs of cluster formation 
and for a wide range of masses, which would 
provide a larger lever arm for comparisons with 
models. They would also provide robust esti- 
mate of total cluster mass using established ro- 
bust, low-scatter observable mass proxies [|53l . 

SZE observations of hot ICM will play an 
increasingly important role in the next decade 
(see white paper by S. Golwala et al.), as a com- 
bination of new sensitive interferometers and 
single-dish telescopes with large, multi-band ~ 
mm focal plane detector arrays should be com- 
ing online. The unprecedented combination of 
high angular resolution (~arcsecond) and spa- 
tial dynamic range of the new planned inter- 
ferometers will allow detailed imaging of the 
ICM thermal pressure within inner regions of 



clusters (6 < 1') with resolution well matched 
to that of the Chandra X-ray observations. 
Such high-resolution observations will also be 
able to estimate the contribution of background 
sources to the SZE signal, thereby eliminating 
one of the main sources of systematic uncer- 
tainty. Continuing development and deploy- 
ment of SZ receivers with large bolometric de- 
tector arrays (> 10 4 detectors) on telescopes 
with wide fields is critical to map out thermal 
pressure of the ICM to the virial radius and be- 
yond for representative samples of clusters. 

The combination of X-ray and new gener- 
ation of SZE observations in the next decade 
would be particularly powerful in mapping the 
distribution and thermal properties of the hot 
ICM out to large radii and in cross-checking 
for systematics. Such observations would de- 
cisively constrain or exclude various scenarios 
of galaxy formation and feedback in clusters, as 
well as plasma effects, such as helium sedimen- 
tation DM EH [561 mi. 

Mapping the cold gas. Mass fractions of cold 
gas are small compared to stars and hot gas, 
but their measurement is important for under- 
standing just how multiphase the ICM really 
is. Such measurements will also provide infor- 
mation on the effectiveness of feedback from 
galaxies and AGN by constraining how much 
of the central gas actually cools. During the 
next decade ALMA will revolutionize measure- 
ments of molecular gas in clusters over a wide 
range of redshifts, while Square Kilometer Ar- 
ray pathfinder missions can be used to map dif- 
fuse HI gas in distant clusters. 

Mapping the total mass distribution. There 
are now several well-established techniques for 
deriving the distribution of total mass in clus- 
ters. Accurate gas distribution and temperature 
profiles obtained with existing and future X-ray 
observatories can provide estimates of the to- 
tal mass profile for hundreds of relaxed clus- 
ters at different redshifts. At the same time, the 
availability of high-resolution SZ observations 
would allow application of such techniques to 



more distant clusters and provide useful cross- 
checks for lower z clusters ll56ll . for which both 
X-ray and SZ data will be available. Additional 
checks will come from measurements of strong 
lensing masses in cluster cores Il58ll59ll , as sta- 
tistical samples of arcs emerge from the upcom- 
ing wide area photometric surveys with follow 
up from space. 

Wide-field deep imaging from large ground- 
based telescopes, but especially by imaging 
cameras on space missions, will allow strong 
and weak lensing measurements of total mass 
profiles (with potential great synergy between 
the two). Such measurements for large samples 
of clusters will provide independent mass mea- 
surements and will allow statistical calibration 
of masses derived from hydrostatic equilibrium 
analyses of X-ray and SZ data. 

3.2 Theoretical modeling. The concerted, 
multi-wavelength observational campaign 
outlined above will need to be accompanied 
by improvements in theoretical modeling of 
cluster baryons. Specific, robust theoretical 
predictions of the properties of cluster galaxies, 
profiles of ICL, and gas density distribution and 
thermal structure of the hot ICM are required. 
Simulations are not yet at the point where 
they can reliably predict properties of cluster 
galaxies or faithfully model various aspects 
of galaxy formation and ICM evolution. The 
focus will therefore be on exploring a wide 
range of plausible scenarios and different 
physical processes: efficiency of star formation 
and its dependence on environment, efficiency 
of stellar and AGN feedback, generation and 
evolution of non-thermal components in the 
form of cosmic rays and magnetic fields (see 
white papers by Myers et al. and Rudnick et 
al.), etc. Progress is particularly needed in 
understanding the role of collisionless effects 
and deviations from equilibrium in plasma. 
This will be essential for reliable interpretation 
of hot gas fraction measurements, especially in 
the cluster outskirts. 



At the same time, advances in available com- 
puting power should allow for larger simula- 
tions which will better resolve small-scale pro- 
cesses during galaxy formation. Note that cos- 
mological simulations of clusters are particu- 
larly difficult in this respect, as scales of kpc 
(or smaller) need to be resolved to model the 
ISM and star formation in galaxies reasonably, 
while modeling volumes in excess of 100 Mpc. 
To simply resolve the ICM within a virial radius 
of a massive cluster uniformly to 1 kpc scale re- 
quires ~ 10 11 resolution elements, well beyond 
capability of the current simulations. Modern 
methods overcome this daunting requirement 
by employing either Lagrangian or adaptive 
mesh refinement techniques, which concentrate 
resolution elements in small fractions of the 
volume in cluster center and around densest re- 
gions of cluster galaxies. Resolving small-scale 
turbulence, shocks, cold fronts, and other po- 
tentially important phenomena over large vol- 
ume will challenge the capabilities of these 
techniques. The wide range of included phys- 
ical processes also significantly increases re- 
quirements for memory, storage, and band- 
width capabilities. This means that substantial 
tera- and peta-scale computing facilities will be 
required for successful work in this direction 
and argues for continuing investment in such fa- 
cilities. Investment into funding computer sci- 
ence research is also needed, as it would help 
ensure that simulation codes can effectively use 
available computing platforms. 
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